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High-field Phase Diagram and Spin Structure of Volborthite Cu 3 V 2 07(OH)2-2H20 
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We report results of V NMR experiments on a high-quality powder sample of volborthite 
Cu3V20t(OH)2-2H20, a spin-1/2 Heisenberg antiferromagnet on a distorted kagome lattice. 
Following the previous experiments in magnetic fields B below 12 T, the NMR measurements 
have been extended to higher fields up to 31 T. In addition to the two already known ordered 
phases (phases I and II), we found a new high- field phase (phase III) above 25 T, at which 
a second magnetization step has been observed. The transition from the paramagnetic phase 
to the antiferromagnetic phase III occurs at 26 K, which is much higher than the transition 
temperatures from the paramagnetic to the lower field phases I (B < 4.5 T) and II (4.5 < B 
< 25 T). At low temperatures, two types of the V sites are observed with different relaxation 
rates and line shapes in phase III as well as in phase II. Our results indicate that both phases 
II and III exhibit a heterogeneous spin state consisting of two spatially alternating Cu spin 
systems, one of which exhibits anomalous spin fluctuations contrasting with the other showing 
a conventional static order. The magnetization of the latter system exhibits a sudden increase 
upon entering into phase III, resulting in the second magnetization step at 26 T. We discuss 
the possible spin structure in phase III. 

KEYWORDS: kagome lattice, frustration, volborthite, Cu3V 2 07(OH)2-2H 2 0, NMR, spin structure, phase 
diagram, high-field 



1. Introduction 

The kagome lattice, a two-dimensional (2D) network of 
corner-sharing equilateral triangles, is known for strong 
geometrical frustration. In particular, the ground state 
of the quantum 5 = 1/2 Heisenberg model with a 
nearest-neighbor antiferromagnetic (AF) interaction on 
the kagome lattice is believed to show no long-range 
magnetic order. Theories have proposed various ground 
states such as spin liquids with no broken symme- 
try with 1 ' or without 2 - 1 a spin-gap or symmetry break- 
ing valence-bond-crystal states. 3 ) Real materials, how- 
ever, have secondary interactions, which may stabilize 
a magnetic order. Effects of the Dzyaloshinsky-Moriya 
(DM) interaction, 4 -' spatially anisotropic exchange inter- 
actions, 5 ~ 7 ) and longer range Heisenberg interactions 8 > 
have been theoretically investigated. In the classical sys- 
tem, it is known that ordcr-by-disorder effect 9 - 1 or DM 
interaction 10 ) favors long-range order with the \/3 x \/3 
pattern or the Q = propagation vector, respectively. 

On the experimental side, intensive efforts have been 
devoted to synthesize materials containing 5=1/2 spins 
on a kagome lattice. n ~ 14 ) Candidate materials known to 
date, however, depart from the ideal kagome model in one 
way or another, such as disorder, structural distortion, 
anisotropy, or longer range interactions. For example, 
herbcrtsmithitc ZnCu3(OH)6Ci2 realizes a structurally 
ideal kagome lattice 12 - 1 and exhibits no magnetic order 
down to 50 mK. 15 ' However, chemical disorder replac- 
ing 10% of Cu 2+ spins with nonmagnetic Zn 2+ ions 16 - 1 
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significantly disturbs the intrinsic properties at low tem- 
peratures. 17, 1S ) 

Volborthite Cu 3 V20 7 (OH)2-2H 2 has distorted 
kagome layers formed by isosceles triangles. Conse- 
quently, it has two inequivalent Cu sites and two kinds 
of exchange interactions (J and J' ) as shown in Fig. 1. 
The Cu2 sites form linear chains, which arc connected 
through the Cul sites. The magnetic susceptibility x 
obeys the Curie- Weiss law x = C/{T + Ow) above 
200 K with Ow = 115 K, exhibits a broad maximum 
at 20 K, and approaches a finite value at the lowest 
temperatures. 11 ' 19 - 1 An unusual magnetic transition 
has been observed near 1 K, 20 ~ 23 ) which is much 
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Fig. 1. (Color online) Schematic structure of volborthite pro- 
jected onto the a-b plane. The H and O sites are not shown. 
The V sites are located both below and above the Cu kagome 
layers. The upper and lower V sites are related by inversion with 
respect to the Cu sites, hence all V sites are equivalent. 
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lower than 9w, consistent with strong frustration in a 
kagomc lattice. However, a recent density functional 
calculation 24 ) proposed that the frustration is attributed 
to the competition between a ferromagnetic J and an 
antiferromagnetic J" between second neighbors along 
the chain (Fig. 1) rather than the geometry of a kagome 
lattice. Thus the appropriate spin model for volbor- 
thite has not been settled yet. Recently, anomalous 
sequential magnetization steps were reported in a high 
quality polycrystalline sample at 4.3, 25.5, and 46 T. 19 ) 
Furthermore, a magnetization plateau was observed at 
2/5 of the saturation magnetization above 60 T, 25 ) in 
spite of the theoretical prediction for a plateau 26 - 1 or a 
ramp 27 ) at the 1/3 of the saturation magnetization in 
an isotropic kagome lattice. 

The magnetic transition at T* ~ 1 K has been de- 
tected by 51 V NMR, 21 ' 22 ) muon spin relaxation, 2 °) and 
heat capacity 23 ^ experiments. The NMR measurements 
revealed a sharp peak in the nuclear relaxation rate 1 /T\ 
and broadening of the NMR spectrum due to develop- 
ment of spontaneous moments for the magnetic field B 
below 4.5 T. 22 ' A kink was also observed in the heat 
capacity. 23 ) Recently, development of short range spin 
correlation was detected by neutron inelastic scattering 
experiments. 28 ) However, the low T phase, which we call 
phase I, shows various anomalies incompatible with a 
conventional magnetic order. 22 ) The NMR line shape 
is not rectangular but can be fit to a Lorentzian, sug- 
gesting a state in which the amplitude of the static mo- 
ment has spatial modulation, such as a spin density wave 
(SDW) state or a spatially disordered state. The behav- 
ior 1/Ti on T provides evidence for dense low-energy ex- 
citations. The spin-echo decay rate I/T2 is anomalously 
large, pointing to unusually slow fluctuations. 

Above 4.5 T, at which the first magnetization step 
occurs, another magnetic phase appears with a different 
spin structure and dynamics. 22 ' 29 ) In this phase, which 
we call phase II, the NMR results show coexistence of two 
types of the V sites with different values of 1 /T2 and line 
shapes. 29 ) The sites with large I/T2 have similar features 
to the V sites in phase I. They are strongly coupled to 
the Cu spins with a spatially modulated structure and 
large temporal fluctuations. On the other hand, the sites 
with small I/T2 show a rectangular spectral shape, which 
is compatible with a conventional order having a fixed 
magnitude of the ordered moment, such as a Neel state 
or a spiral order. The numbers of these two sites are 
nearly equal in a wide range of B and T. Based on these 
observations, a heterogeneous spin state, in which two 
types of Cu moments form a periodic structure, has been 
proposed. 29 ) 

In this paper, we report results of the NMR experi- 
ments in higher magnetic fields up to 31 T. The phase 
diagram determined from the previous 22,29 ) and the 
present experiments is shown in Fig. 2. A new high- 
field phase defined as phase III is found above 25 T, at 
which the second magnetization step occurs. The transi- 
tion from the paramagnetic phase to the phase III takes 
place at 26 K, which is much higher than the transition 
temperatures from the paramagnetic to phases I and II. 
At low temperatures, two types of the V sites are ob- 
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Fig. 2. The phase diagram of volborthitc proposed from the pre- 
vious 19 ' 22 ' and present studies. The squares (triangles) represent 
the phase boundaries determined from the variation of the NMR 
spectra as a function of temperature (magnetic field). The cir- 
cles and stars indicate the peak of 1/Ti 22 ' and magnetization 
steps, 19 ' respectively. The lines are guide to the eyes. There is a 
finite range of field, where phase II and III coexist. 



served in phase III in a similar way as in phase II, indi- 
cating that the high-field phase III exhibits also a hetero- 
geneous spin state. Our results show that those Cu spins 
which arc dominantly coupled to the V sites with small 
I/T2 are responsible for the increase of magnetization at 
the second step. We propose that the second magneti- 
zation step is ascribed to the ferromagnetic alignment 
of these Cu spins and discuss possible spin structure in 
phase III, with particular reference to the theories on the 
anisotropic kagome lattice in the limit J 3> J'. 5,6 ) 

2. Experiment 

The 51 V NMR measurements have been performed on 
a high-quality powder sample prepared by the method 
described in ref. 19. The data at high magnetic fields 
above 18 T were obtained by using a 20 MW resistive 
magnet at LNCMI Grenoble. The NMR spectra were 
recorded at fixed resonance frequency (vq) by sweeping 
magnetic field (B) in equidistant steps and summing the 
Fourier transforms of recorded spin-echos, which were ob- 
tained by using the pulse sequence 7r/2 — r — 7r/2. 30 ) They 
are plotted against the internal field B lnt = uq/'J—B. The 
NMR spectra are labeled by the reference magnetic field 
Bq = ^0/7- We determined 1/Ti by fitting the spin-echo 
intensity M(t) as a function of the time t after a comb 
of several saturating pulses to the stretched exponential 
function M(t) = M eq - A/ cxp{-(t/Ti)' 3 }, where M eq is 
the intensity at the thermal equilibrium. This functional 
form was used to quantify the inhomogeneous distribu- 
tion of 1/Ti. When the relaxation rate is homogeneous, 
the value of (3 is close to one. 
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Fig. 3. (Color online) Temperature dependence of the NMR spec- 
tra at Bo = 30.1 T. The inset shows an example of calculated 
powdcrpattcrn broadened by a non-axial anisotropic magnetic 
shift. 



3. Results and Discussion 

3.1 Temperature dependence of NMR spectra in high 
magnetic fields 

Figure 3 shows the temperature dependence of the 
51 V NMR spectra at B = 30.1 T. The spectra above 
30 K are compatible with a powder pattern due to 
anisotropic magnetic shifts (inset of Fig. 3). The width 
due to the quadrupolar interactions is about 0.02 T, 22 ) 
much smaller than the observed width in this field range. 
As temperature decreases, the spectral shape begins to 
change near 26 K. Below 20 K, the spectra show two 
peaks. Such a structure, which cannot be explained by 
the anisotropic paramagnetic shift, indicates an inter- 
nal field due to spontaneous antiferromagnetic (AF) mo- 
ments. Therefore, our result indicates an AF transition 
near 26 K at 30 T. Similar measurements at 25.6 T in- 
dicate the transition near 22 K. As shown in Fig. 3, the 
spectral width rapidly increases with decreasing temper- 
ature below 15 K. The spectral shape becomes insensitive 
to the temperature below 4 K, where three broad peaks 
labeled as A, B, and C are observed. 

The transition temperature of 26 K is much higher 
than those between the paramagnetic phase and phase I 
or II below 12 T. 22 ) Therefore, this suggests a new mag- 
netic phase at high fields which is distinct from phases I 
and II. Indeed, the magnetization data shows the second 
step near 26 T, indicating a field-induced magnetic tran- 
sition. This new high-field phase will be called phase III, 
as shown in Fig. 2. 

Figure 4 shows the temperature dependences of 1/Ti 




Fig. 4. (Color online) Temperature dependences of 1/Ti and the 
strech exponent /3 at 30.1 T. 
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Fig. 5. (Color online) Recovery curves M(t) measured at posi- 
tions A, B, and C in Fig. 3 at 1.3 K. The solid lines show 
the fitting to the stretched exponential function M(t) = M eq — 
M exp{-(t/Ti)' 3 }. 



and the strech exponent (3 at 30.1 T. In the tempera- 
ture range between 15 and 50 K, 1/Ti is independent of 
temperature. In spite of the clear change in the spectral 
shape, there is no significant anomaly in 1/Ti near 26 K. 
Let us recall that while the transition from the param- 
agnetic phase to phase I is marked by a clear peak of 
1 jT\ . only a small anomaly in 1 jT\ was observed at the 
transition from the paramagnetic phase to phase II. 22,29 ' 
Thus the dynamic signature at the magnetic transition 
becomes gradually obscured for higher-field phases. How- 
ever, it is still puzzling that, unlike in phase II, 1/Xi no 
longer decreases below the magnetic ordering tempera- 
ture in phase III. 

As shown in Fig. 4, f3 decreases with decreasing tem- 
perature below 10 K, indicating more inhomogeneous 
distributions in 1/Ti. With further decreasing temper- 
ature, the distribution of 1/Ti becomes too large to de- 
termine the representative value of 1/Ti. In addition, the 
relaxation rates depend on the spectral position at low 
temperatures. Figure 5 shows the recovery curves M{t) 
measured at positions A, B, and C in Fig. 3 at 1.3 K. 
The intensity at A is recovered within 0.1 s, while the in- 
tensity at C does not saturate even at 10 s. These curves 
can be fit to the stretched exponential function, and we 
obtain 1/Ti = 311 s" 1 and /3 = 0.49 for A and 1/Ti 
= 2.9 s _1 and /3 = 0.33 for C. The recovery curve at 
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Fig. 6. (Color online) Magnetic field dependence of the NMR 
spectra at 0.4 K. Straight lines arc guide to the eye to follow 
different components (peaks) of the spectra. 



B shows an intermediate behavior between the curves 
measured at A and C. Consequently, M(t) at B does 
not fit to the stretched exponential function. These re- 
sults show that the spectra at low temperatures consist 
at least of two components. This is quite similar to the 
behavior observed in phase II, where the previous NMR 
measurements show the coexistence of two components 
characterized by the large and small values of I/T2. 29 ' 
The two component behavior will be discussed in detail 
in section 3.3. The onset of rapid spectral broadening 
(Fig. 3) and the strongly inhomogeneous distribution of 
1/Ti (Fig. 4) indicate that there may be another tran- 
sition or crossover at 10 - 15 K. This is shown by the 
dotted line in Fig. 2. 

3.2 Magnetic field dependence of NMR spectra 

Figure 6 shows the field dependence of the 51 V NMR 
spectra at 0.4 K. The center of gravity of the spectra 
shifts to larger values of B lnt with increasing magnetic 
field. This shift corresponds to the increase of the mag- 
netization. On the other hand, the width and shape of 
the spectrum are related to the magnetic structure. As 
shown in Fig 6, the spectral shapes are almost indepen- 
dent of the magnetic field in the ranges of 18 - 21 T and 
above 28 T. However, the spectra above 28 T have a dif- 
ferent shape from the spectra below 21 T. This indicates 
that the transition from phase II to phase III occurs in 
the magnetic field region between 21 and 28 T, and the 
magnetic structure in phase III is significantly different 
from that in phase II. In Fig. 6, we observe that the 
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Fig. 7. (Color online) Magnetic field dependence of the center of 
gravity of the spectrum M\. The experimental temperature is 
0.3 K for the data below 15 T and 0.4 K for the data above 
15 T. The open and solid circles represent Mi for the slow and 
fast components, respectively (see the text). 



broad peak indicated by the dashed line is present in all 
the spectra, shifting slightly to higher B lnt with increas- 
ing magnetic field. In phase II, an additional peak and 
a shoulder indicated by the dotted lines are observed at 
lower and higher B lnt , respectively. They gradually dis- 
appears above 22 T. Instead, two additional peaks indi- 
cated by the solid lines are observed in phase III. These 
two peaks gradually develop above 24 T. As the magnetic 
field increases from 22 to 28 T, there is a gradual transfer 
of the spectral weight of the peaks denoted by the dotted 
lines to those denoted by the solid lines, indicating the 
coexistence of phases II and III. 

In order to investigate the phase boundary in more 
detail, we examine the first and second moments of the 
spectra defined as, 



Mi 



Mo 



BintI(Bint)dBi, 



(Bint - Mi) 2 I{B int )dB int 



(1) 



where I(Bi nt ) is the NMR spectrum normalized as 
J I(Bi nt )dBi nt = 1. Figure 7 shows the magnetic field 
dependence of Mi. The data below (above) 15 T were 
taken at 0.3 K (0.4 K). At these sufficiently low tem- 
peratures, difference of 0.1 K makes no change in the 
spectral shape. As shown in Fig. 7, the data below 22 T 
lie on a straight line through the origin. Above 22 T, Mi 
departs from the straight line. The magnetic field de- 
pendence of Mi in Fig. 7 reproduces the magnetization 
curve 19 - 1 including the second magnetization step at 26 
T. 

Figure 8 shows the magnetic field dependence of the 
line width, i.e. square root of the second moment M2, 
at various temperatures. At 0.3 K, \fM-2 decreases with 
increasing B below 4.5 T (phase I). For the range of mag- 
netic field 5 - 22 T and at 0.3 - 0.4 K, is almost 
independent of B, indicating that magnetic structure re- 
mains unchanged over the entire field range of phase II. 
Above 28 T, \pKl2~ becomes independent of the magnetic 
field again. Thus the region above 28 T should belong to 



J. Phys. Soc. Jpn. 



Full Paper 



Author Name 5 



0.10 





A* 




1 T 


▼ 


o 


6 


■ 


■ 


8.5 




A 


11.5 
14 


▼ 
■ 


▼ 


16 


* 


□ 


18.3 



▼ . 

° cP B 



o 



n 



2 3 

J(K) 



a 

B 



Fig. 8. (Color online) Magnetic field dependence of square root 
of the second moment M2. The open and solid circles represent 
the data at 0.3 and 0.4 K, respectively. The solid triangles and 
squares represent the data at 1.35 and 4.2 K, respectively. 



Fig. 9. (Color online) Temperature dependences of FWHM at 
various magnetic fields. The data below 11.5 T are from in ref. 
22. 



phase III. For the range of B between 22 and 28 T, 
increases monotonically with increasing magnetic field. 
This field region is likely to correspond to the broadened 
transition due to anisotropy in the powder sample. The 
similar broadened transition between phases I and II is 
also observed around 4.5 T with the width of the coexis- 
tence region of about 1 T. 22 ' In both cases, the width of 
the transition is about 20 - 25% of the center field of the 
transition. A likely origin of the broadened transition is 
the anisotropy of g factor, which is estimated to be about 
17% from the electron spin resonance measurements. 31 ' 

As shown in Fig. 8, \fM% at 0.4 and 1.35 K show sim- 
ilar behavior above 18 T, indicating that the transition 
field is unchanged at these temperatures. On the other 
hand, \JM^ at 4.2 K lie on a straight line through the 
origin below 18 T. This behavior indicates a paramag- 
netic state, where the internal field is induced by the 
external magnetic field. However, y ' M2 departs from the 
straight line above 22 T, indicating the appearance of 
a spontaneous internal field. Above 28 T, \JM-i becomes 
independent of the magnetic field. Therefore, for the field 
above 28 T, there must be an AF order even at 4.2 K. 
Thus obtained phase boundaries arc plotted by the solid 
triangles in Fig. 2. 

In order to establish the phase boundary between 
phase II and the paramagnetic phase completely, the 
temperature dependence of the full width at half maxi- 
mum (FWHM) of the spectra was measured as shown in 
Fig. 9 at various magnetic fields. Here, the data below 
11.5 T are the same as those presented in ref. 22. The 
transition temperatures are determined by the onset of 
line broadening. Above 8.5 T, the onset is insensitive to 
the magnetic field as indicated by the arrow in Fig. 9. 
The transitions observed by our NMR measurements are 
summarized in Fig. 2. Phases I and II are characterized 
by a weak field dependences of the transition tempera- 
ture. Phase III is characterized by much higher transition 
temperatures compared with phases I and II. 

3.3 Spin structure 

To consider spin structures in the low-temperature 
phases of volborthitc, we begin with general discussion 



on the powder-pattern NMR spectra for various mag- 
netic states. Figure 10 shows the schematic illustration 
of the typical powder-pattern for three different antifer- 
romagnctic states. The origin of the horizontal axis is set 
at the center of gravity of the spectrum. Fig. 10(a) shows 
a Gaussian spectrum which is expected if the magnitude 
of the internal field is randomly distributed. Fig. 10(b) 
shows a spectrum for a spin density wave (SDW) state, 32 ) 
where the magnitude of the internal field has a sinusoidal 
modulation. Fig. 10(c) shows a rectangular spectrum, 
which is expected when the internal field has a fixed mag- 
nitude. Examples of this case include a simple Neel state 
with a two-sublattice structure and a spiral order with an 
isotropic hyperfine coupling. For both cases (b) and (c), 
we assume that the direction of the internal field is ran- 
domly distributed with respect to the external field, i.e. 
spin-flop transition does not occur. As shown in Fig. 10, 
both the SDW and Gaussian distribution give a central 
peak. They differ only in the way the intensity vanishes 
on both sides. Since NMR spectra of powder samples are 
determined by the distribution of the magnitude of B mt , 
it may be difficult to distinguish from the NMR spectra 
alone between an SDW order, where i?i n t has a coherent 
spatial dependence as cos(Q • r), and a random distri- 
bution of Bi nt , in particular if the distribution is not a 
Gaussian. On the other hand, it is easy to distinguish 
those ordered states with spatial modulation in the mag- 
nitude of moments such as SDW from the ones with a 
fixed magnitude of moments. The spectra observed in 
phase I belongs to the former case as reported in ref. 22, 
i.e. the magnitude of the internal field is spatially mod- 
ulated. 

In phase II, the previous NMR results show coexis- 
tence of two types of V sites, V/ and V s sites, which 
are characterized by the large (fast) and small (slow) re- 
laxation rates, respectively. 29 ) The spin-echo decay curve 
in phase II was well fit to the two-component function 
J(t) = A / cxp(-2r/T 2/ ) + A s cxp(-2r/T 2s ). The coeffi- 
cients A f and A s are proportional to the number of V/ 
and V s sites within the frequency window covered by 
the exciting rf-pulse. Figure 11(a) shows the spectra for 
the V/ and V s sites at 6 T and 0.3 K reported in ref. 
29. These spectra are obtained by plotting Af and A s at 
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Fig. 11. (Color online) (a) NMR spectra for the fast (solid circles) 
and slow (open circles) components at 0.3 K in the magnetic 
field of 6.0 T. (b) NMR spectra at r = 10 (thick solid line) and 
110 fis (thick dashed line) in the magnetic field of 30.1 T. The 
red (blue) dotted line represents the NMR spectra for the fast 
(slow) component. 



Fig. 10. (Color online) Typical powder-pattern spectra for three 
cases, (a) Gaussian distribution in the magnitude of the internal 
field , (b) a SDW state where the internal field has a sinusoidal 
modulation, and (c) an Neel state with a simple two-sublatticc 
structure. 



various frequencies v$ with the fixed B = 6.0 T, after the 
two-component fitting has been done. The spectrum of 
the V/ sites resembles a Gaussian. Therefore, the V/ sites 
in phase II have a similar feature to the V sites in phase 
I. These V sites should be surrounded by Cu moments 
with spatially varying magnitude. On the other hand, the 
spectrum of the V s sites shows a rectangular-like shape, 
indicating that the internal field at the V s sites is more 
homogeneous, compared with that at the V/ sites. The 
spin structure of the V s sites can be a two-sublatticc 
Neel order or a spiral order with an isotropic hypcrfinc 
coupling. When the hypcrfinc coupling tensor is slightly 
anisotropic, a spiral order leads to a small modulation of 
the internal field, resulting in a trapezoidal spectrum. In 
reality, this may be difficult to distinguish from a rectan- 
gular one. Thus we cannot select a unique spin structure. 

The multi-component behavior is observed also in 
phase III by the 1/Ti measurements at low temperatures 
as shown in Fig. 5. These components also exhibit differ- 
ent spin-echo decay rates. In Fig. 11(b), two spectra at 
30.1 T and 2.5 K are displayed, one for r = 10 //s and 
the other for r = 110 /is. The spectrum for r = 10 /is 
is the same as the one shown in Fig. 3, and the peaks 



are labeled as A, B, and C. The intensity of the peak 
A markedly decreases with increasing r from 10 to 110 
/its. That is, the A component has large 1/T 2 as well as 
large 1/Ti. The left half of the peak A of the spectrum 
for t = 10 /xs can be fit to a Gaussian labeled as "fast" , 
which disappears in the spectrum for r = 110 /is. By sub- 
tracting the Gaussian component from the spectrum for 
t = 10 /us, we obtain a rectangular component shown by 
the dotted line labeled as "slow" . Then, the "slow" com- 
ponent is quite similar to the spectrum for r = 110 /is, 
indicating that the two-component description is valid 
also in phase III. 

The common feature of phases II and III are two types 
of V sites, the V/ and V s sites, of similar characteristics. 
At the V/ sites, the spectral shapes resembles a Gaussian 
and the relaxation rates are large. On the other hand, the 
V s sites have a rectangular-like spectral shape with rela- 
tively small relaxation rates. However, the average values 
of Bi nt for the two V sites are quite different for phases 
II and III as shown in Figs. 11(a) and (b). In going from 
phase II to phase III, the center of gravity of the whole 
spectrum shifts to higher values of B lnt due to increase of 
the magnetization. What is more interesting is that the 
relative positions of the two components change substan- 
tially between phases II and III. In phase II, the centers 
of gravity of the fast and slow components are located 
almost at the same position, as shown in Fig. 11(a). This 
shows that the two types of Cu spins, each of which is 
dominantly coupled to either Yf or V s sites, have al- 
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most the same field-induced averaged magnetization. On 
the other hand, in phase III, the center of gravity of the 
slow component is located at a much larger value of _Bi, lt 
than that of the fast component. Therefore, the averaged 
magnetization coupled to the V s sites is much larger than 
that coupled to the V/ sites in phase III. The centers of 
gravity of the fast and slow components above 28 T are 
shown in Fig. 7. The former is on the straight line extrap- 
olated from phase II, while the latter is located at much 
larger i?i n t- Therefore, we conclude that the Cu spins 
coupled to the V s sites are responsible for the increase 
of magnetization at the second step. 

The center of gravity Mi is related to the magneti- 
zation M by the relation M = Mi /Am, where Am is 
the coupling constant determined in the paramagnetic 
phase Am = 0.77 T / )j,b 21 ' 22 ' Indeed, the center of grav- 
ity of the whole spectrum at 31 T, Mi = 0.165 T, gives 
the magnetization of 0.21 \xb in good agreement with 
the measured magnetization of 0.22 hb- 19 ^ By using this 
relation, the averaged magnetization coupled to the Vj 
(V s ) sites, Mf ast (M s i ow ), at 31 T arc estimated to be 
0.16 (0.32 hb)- The integrated intensities of the spectra 
of the fast and slow components are almost the same. 
This means that numbers of the two sites Nf and N s 
are nearly equal. More precise determination of the ratio 
Nf : N s would require precise knowledge of the frequency 
dependence of 1/Ti and 1/T 2 . 

In phase II, a heterogeneous spin state was pro- 
posed, 29 - > in which the V/ and V s sites spatially alter- 
nate, because of magnetic superstructure on Cu sites. 
Arguments were given in ref. 29 why this is more likely 
than the macroscopic phase separation. First, the ratio 
Nf : N s should vary as a function of B and T, if the 
V/ and V s sites belong to distinct phases. Experimen- 
tally, however, Nf and N s are nearly equal irrespective 
of B and T. Second, the V/ and V s sites show similar 
T-dependence of 1/T 2 . 29) This means that the two types 
of Cu spins coupled to the V/ and V s sites, respectively, 
are interacting in a certain way to establish a common 
T-dcpendcncc. Again, this is not possible if they belong 
to different phases. 

Although it is difficult to determine the superstructure 
from the NMR data on a powder sample alone, plausible 
structures were discussed in the previous paper as fol- 
lows. 29 -* The most simple way to generate two types of V 
sites with the same abundance is to form a superstructure 
doubling the size of the unit cell. On the ideal kagome 
lattice, the three directions ki, k 2 , and ki + k 2 (Fig. 12 
a) are all equivalent. However, in volborthite, the struc- 
ture is uniform along k 2 and inequivalent Cul and Cu2 
sites alternate along ki. Then it is more likely that the 
superstructure develops along ki. Consequently, each of 
the Cul and Cu2 sites are divided into two inequivalent 
sites Culp, Cul CT and Cu2 K , Cu2 A , as shown in Fig. 12(a). 
It follows that there are two types of the V sites, V Q and 
V/3, as shown in Fig. 12(a). They should correspond to 
the V/ and V s sites. It should be emphasized that dis- 
tinct spin states on the crystallographically inequivalent 
Cul and Cu2 sites do not account for our observation of 
two different V sites that are otherwise crystallographi- 
cally equivalent. Our results can be accounted for only 



by a symmetry breaking superstructure. 

Since a V site is located approximately above or below 
the center of Cu hexagon of the kagome lattice, the dom- 
inant source of the hypcrfinc field at V nuclei should be 
confined within the six Cu spins on a hexagon. Because of 
the distorted structure of volborthite, there are three dis- 
tinct hyperfine coupling tensors A , A&, and A c as shown 
in Fig. 12(a). The coupling tensors to the other three 
spins A^, A! b , and A' c are generated by the mirror reflec- 
tion perpendicular to the b axis at the V site. A and A' 
are the same hyperfine tensors but oriented differently, 
and they will be effectively different only if the tensor 
is anisotropic. The coupling constant Am determined in 
the paramagnetic state should be equal to the average of 
the diagonal components of the sum of these six tensors. 
The hyperfine fields at the V Q and Yp sites are written as 
B a = (A s K +A / s^)+(A 6 s p +A^sp+(A c SA+A' c s' A ) and 
= (A a s A +A^s' A ) + (A 6 s (T +A , b s;) + (A c s K +A^^), re- 
spectively, where s £ and s' e denote two neighboring spins 
on the same type of sites (i.e., e stands for Cul p , Cul,j, 
Cu2 K , or Cu2 A ). 

Both 1 jT\ and 1 /T 2 are determined by the time corre- 
lation function of the hyperfine field. 33 ) Since the hyper- 
fine coupling to the uniform magnetization Am is largely 
isotropic, 22 ' we expect this is also the case for the individ- 
ual coupling tensors A a , Aft, and A c . If they have similar 
magnitudes, the coupling to the Cu2 sites is identical for 
the two V sites, and the different relaxation rates must 
be ascribed to the difference in the time correlation func- 
tions of the Cul spins s p and s a . The model proposed 
in ref. 29 is shown in Fig. 12(b). The Culo- sites develop 
a long range magnetic order with a uniform magnitude 
of moments. The spin fluctuations at these sites should 
have a small amplitude, leading to the slow relaxation 
at the V s sites. On the other hand, the Cul p sites show 
a modulation in the magnitude of the ordered moments. 
The absence of fully developed static moments would al- 
low unusually slow fluctuations responsible for the large 
1/T 2 at the V/ sites. 

In phase II, the spectra of the V/ and V s sites have 
almost the same Mi, which is proportional to the field. 
Therefore, all Cu sites have both an AF component and 
a field- induced uniform component. In phase III, Mi is 
largely shifted only at the V s sites. This shift can be 
explained by ferromagnetic alignment of the Culo- sites 
as shown in the right panel of Fig. 12(b). As shown in 
Fig. 11, the width of the spectrum of the V s sites in phase 
III is narrower than that in phase II. This decrease of the 
width is also understood by the ferromagnetic saturation 
of the Cul CT sites. Because the ferromagnetic moment is 
always parallel to the magnetic field, the spectral width is 
determined by the anisotropy of the hyperfine coupling. 
Indeed, the spectral shape of the V s sites in phase III is 
similar to a paramagnetic powder pattern. In contrast, 
the width in phase II should be dominantly determined 
by the AF component. On the other hand, the position 
and the width of the spectrum of the V/ sites do not 
change substantially at the transition between phases II 
and III. This indicates that the Cul p spins keep spatial 
modulation of the magnitude of moment in phase III. We 
should stress that the width of the V s spectrum is much 
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Fig. 12. (Color online) Possible magnetic structures in volbor- 
thitc projected onto the a-b plane, (a) Superstructure along the 
ki direction provides two incquivalent V sites, V a and Vp, and 
four incquivalent Cu sites, Cul p , Cul CT , Cu2 K , and Cu2^. (b) If 
the coupling tensors A a , Aj,, and A c are largely isotropic and 
have similar magnitudes, different relaxation behavior at the 
and V s sites should be ascribed to the distinct spin states at the 
Cul/, and Cul CT sites, (c) If A a is the dominant coupling, differ- 
ent fluctuations at the Cu2 K and Cu2^ sites lead to distinction 
between V/ and V s . 



smaller than the shift Mi in phase III, consistent with the 
largely isotropic hyperfine coupling in the paramagnetic 
phase. 22 ) 

Our results show interesting correspondence with the 
theory on the anisotropic kagome model in the limit 
J ^ J'. 5,6 ) Schnyder, Starykh, and Balents studied the 
ground state of this model in zero magnetic field. They 
found that the Cul sites develop cither a ferromagnetic 
order or a spiral order with a long wave length with the 
full moment, while the Cu2 sites have only a small or- 
dered moment. 5 ^ The spin fluctuations at the Cu2 sites 
associated with the large energy scale ~ J should also 
lead to a small contributions to 1/Ti and I/T2 at the 
V sites. The irrelevance of the Cu2 sites for both the 
static internal field and the dynamic relaxation rates at 
the V sites is consistent with our model that the differ- 
ence between the V/ and V s sites are ascribed to the two 
types of the Cul sites. It should be noted, however, that 
the theory does not account for the superstructure which 
produces the different V/ and V s sites. 



The similarity between our model shown in Fig. 12(b) 
and the theoretical prediction on the anisotropic kagome 
model becomes more apparent in high magnetic fields. A 
moderate magnetic field should easily drive the small-Q 
spiral order of the Cul sites into a ferromagnetic state. 
This situation is consistent with our observation of in- 
duced ferromagnetism of Culer spins in phase III. How- 
ever, our results indicate that only a half of the Cul 
sites undergoes ferromagnetic saturation across the sec- 
ond magnetization step at 26 T. It appears most likely 
that the other half of the Cul sites (the Cul p sites) get 
fully polarized across the third step at 46 T. Stoudcn- 
mire and Balents investigated the spin structure of the 
Cu2 sites for the same anisotropic kagome model in high 
magnetic fields, where Cul moments are fully aligned by 
the field. 6 ) They found that the Cu2 sites show a quan- 
tum phase transition from a ferrimagnetic order along 
the field direction to an antiferromagnetic order perpen- 
dicular to the field with increasing field. A remarkable 
aspect of the theory is the prediction for a magnetiza- 
tion plateau at 2/5 of the saturation magnetization in 
the ferrimagnetic region immediately below the transi- 
tion field. 

The 2/5 plateau has been actually observed in the re- 
cent magnetization measurements above 60 T. 25 ) There- 
fore, we further examine if such a scenario is compatible 
with our experimental results. When the Cul sites ex- 
hibit the full moment, (s z ) = 1/2, in the plateau region, 
the Cu2 sites should have the averaged moment (s z ) = 
1/20 in order to account for the total magnetization of 
2/5 of the saturation. If we assume that all the hyperfine 
coupling tensors A a , Af,, and A c are isotropic and have 
the same magnitude, A a = Ai, = A c = A n f/6 = 0.13 
T/hb, the contribution from the four Cu2 sites to B{ nt 
is 0.055 T at all V sites in the plateau region. Here, we 
used the averaged g value of 2.15. 31 ' 34 * 1 This contribu- 
tion to Bi n t will be reduced to 0.028 T at 30 T, assuming 
that the magnetization at the Cu2 sites is proportional 
to B. In addition, the full moments at the Cul^ sites 
produce B- lnt = 0.28 T at the V s sites. The sum of these 
gives i?i n t = 0.30 T, which is reasonably close to the 
observed Mi (= 0.24 T) at the V s sites at 30 T (Fig. 
7). On the other hand, M x = 0.12 T at the V/ sites at 
30 T corresponds to (s z ) — 0.17 at the Cul p sites. The 
summation of the magnetizations at the Cu2, Cul CT , and 
Culp sites amounts to the total magnetization of 0.28 
/is at 30 T, which is in reasonable agreement with the 
experimental magnetization of 0.22 [Ib- Therefore, this 
model is semiquantitatively compatible with the NMR 
and magnetization measurements. 

Let us remark that the V s and V/ sites show ex- 
tremely different relaxation rates in phase III. As shown 
in Fig. 5, the V s sites (line C) show two orders of magni- 
tude smaller 1/Ti than the V/ sites (line A). The spectra 
in Fig. 11(b) for different values of r indicate orders of 
magnitude difference in I/T2. In contrast, the difference 
in I/T2 between the V/ and V s sites is only a factor of 
four in phase II. 29 ' In fact the extremely small relaxation 
rates at the V s sites in phase III is what should be ex- 
pected if the moments at the Cul CT sites are completely 
saturated by magnetic fields and the contribution from 
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the Cu2 sites is heavily suppressed for J ^> J'. On the 
other hand, in phase II, the spins at the Cul CT sites have 
finite fluctuations even though they show a conventional 
magnetic order. The anomalous fluctuations of the Cul p 
spins can then couple to the fluctuations of the Culer 
spins, leading to the similar temperature dependence of 
1/T 2 at the V/ and V s sites. 29 ) 

The above discussion is based on the assumption that 
all the hypcrfinc coupling tensors are largely isotropic 
and have similar magnitudes. However, the low symme- 
try of the volborthite structure may result in large dif- 
ference in the hypcrfinc coupling. In particular, A a and 
A c involve significantly different V-O-Cu hybridization 
paths. If one of them, say A a , is the dominant coupling, 
the different relaxation rates for the two V sites must be 
ascribed to the different dynamics of s>, and s K . An ex- 
ample of this case is shown in Fig. 12(c). However, there 
has been no microscopic theory which supports such a 
scenario. 

4. Summary 

We have presented the 51 V NMR results of the S = 
1/2 distorted kagome lattice volborthite. Following the 
previous experiments in magnetic fields below 12 T, the 
NMR measurements have been extended to higher fields 
up to 31 T. In addition to the already known two ordered 
phases (phases I and II) , we found a new high-field phase 
above 25 T, at which the second magnetization step oc- 
curs. This high-field phase (phase III) has the transition 
temperature of about 26 K, which is much higher than 
those of phase I (1 K) and phase II (1 - 2 K). At low 
temperatures, two types of the V sites are observed with 
different relaxation rates and line shapes in phases II 
and III. Our results indicate that both phases exhibit a 
heterogeneous spin state, where an ordered state with 
non- uniform magnitude of moment with anomalous fluc- 
tuations alternates with a more conventional order with a 
fixed magnitude of the ordered moment. The latter group 
of spins are responsible for the increase of magnetization 
at the second step between the phase II and phase III. 
We proposed a possible spin structure in phase III com- 
patible with the NMR and magnetization measurements. 
Our model has common features with the theories on 
the anisotropic kagome lattice in the limit J 3> J', 5 ' 6 -* 
although the formation of heterogeneous superstructure 
was not predicted by the theories. 
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